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CHAPTER 1 
INTRODUCTION 

Although the development of radar did not mature until World War II, the study of 
electromagnetic scattering, the basis of radar, is almost as old as Maxwell's equations 
themselves. Hertz, experimentally testing the theories of Maxwell in 1886, demonstrated that 
radio waves could be reflected by metallic and dielectric bodies [1]. Sommerfeld's solution 
for the half-plane scattering problem [2] appeared in 1896, and Mie's classic paper on the 
scattering from spheres [3] was published in 1908. In a speech delivered to the Institute of 
Radio Hngineers in 1922 [4j, Marconi said: 

As was shown by Hertz, the electric waves can be completely reflected by conducting 
bodies. In some of my tests I have noticed the effects of reflection and detection of these waves by 
metallic objects miles away. 

It seems to me that it should be possible to design apparatus by means of which a ship 
could radiate or project a divergent beam of these rays in any desired direction, which rays, if 
coming across a metallic object, such as another steamer or ship, would be reflected back to a 
receiver screened from the local transmitter on the sending ship, and thereby, immediately reveal 
the presence and bearing of the other ship in fog or thick weather. 

Today, radar is used not only in ship detection applications as Marconi had predicted, but 
also in air traffic control, remote sensing, satellite detection and tracking, and law enforcement. 
Additionally, radar has a wide variety of military applications: for surveillance and detection of 
friendly and hostile planes, ships, tanks, and missiles; for navigation; and for control and 
guidance of weapons. 

With the increasing effectiveness and sophistication of today's military radar systems, the 
ability to reduce the radar visibility of potential targets has been of importance. There are four 
basic techniques for reducing the radar cross section (RCS) of a target [5]: 
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(i) shaping 

(ii) radar absorbing (or lossy) materials 

(iii) passive cancellation 

(iv) active cancellation. 

Of the four methods, shaping and the use of absorbing materials are the most effective. For 
example, a police radar capable of detecting the large planar surfaces of an Allied moving van 
(Fig. 1.1(a)) 1.5 miles away could not detect the sleek curved surfaces of a Corvette (Fig. 
1.1(b)) 600 feet away [6]. 

One problem of interest in RCS reduction studies is that of the open-ended cavity [7] - 
[12]. Coating the interior of the cavity with lossy materials and shaping the cavities have been 
used in reducing the RCS of cavity structures. This thesis investigates both issues in the study 
of RCS reduction of arbitrary cavity structures. 

In Chapter 2, the effects of interior coating are studied. An open-ended cavity is often 
modelled as a section of waveguide, open on one end and terminated by a flat conducting plate 
on the other. A layer or layers of material lining the waveguide wall serve to alter the modal 
fields in the guide in order to achieve either more or less attenuation for certain modes. With 
the proper choice of coating material and thickness, a significant RCS reduction can be 
achieved using only a single thin layer of coating. However, this reduction is only effective 
over a narrow frequency band at low frequencies (aA ~ 1, where a = radius of the cylinder; 
X = the free-space wavelength). A generalized method is presented in this chapter for 
computing the propagation/attenuation constants for the normal modes of a circular waveguide 
lined with multiple layers of coating. This method is then applied to multilayered coated 
structures to show that a greater effectiveness in the attenuation of the normal modes can be 
achieved over a broadened frequency range. 



3 


ORiGi;^-/;.L PAQE ?S 

OF POOR QUALITY 



Fig. 1 . 1 Comparison of body shapes for (a) Allied moving van, with flat planar surfaces 
and (b) Chevrolet Corvette, with smooth, sleek curves. 
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In Chapter 3, the problem of calculating the RCS of an arbitrary cavity is addressed. 
Traditional methods of analysis for the open-ended cavity include such techniques as the modal 
analysis used in Chapter 2, finite element and difference methods, and moment methods. 
However, the aforementioned techniques become impossible if (i) the shape of the cavity is not 
a perfect sphere or cylinder, (ii) the electrical dimension of the cavity is large, or (iii) the space 
inside the cavity is not homogeneous. These restrictions prevent the realistic modelling of 
physical problems and hinder the study of the effects of shaping on RCS. In this chapter, a 
different strategy for analyzing the open cavity problem, entitled "shooting and bouncing rays" 
(SBR), is presented. A dense grid of geometric optics (GO) rays representing an incident 
plane wave is 'shot' into the cavity through the front aperture and followed as the rays bounce 
from conductors, penetrate through materials, and eventually return to the opening of the 
cavity. An innovative scheme is then used to integrate the aperture field to obtain the scattered 
field. This is the first attempt of its kind at solving partially open structures using the ray 
technique found in the open literature. The SBR method places no restriction on the shape of 
the cavity and permits the RCS computation for any arbitrarily shaped cavity. 

Chapter 4 studies the effectiveness of shaping as an RCS reduction method. The SBR 
method is used to study the RCS of uniform and nonuniform cylindrical cavities. For a large 
(relative to wavelength) flat-plate terminated cylinder, there is a large RCS value at near axial 
incidence due to the reflection from the termination . The near axial RCS of such structures can 
be reduced by a smooth gradual bending of the cylindrical axis. This chapter shows the 
reduction of RCS achieved by applying a longitudinal S-bend to a circular waveguide and a 
triangle-to-circle transition waveguide. The SBR method is currently the only available viable 
method for computing the RCS for such geometries. 

The range of validity of the SBR method is investigated in Chapter 5. The limitations of 
SBR are studied by considering the problem of a plane wave impinging on a semi-infinite 
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parallel-plate waveguide. Two alternative descriptions for the fields propagating in the 
waveguide are the modal analysis and GO ray analysis, used by the SBR method. A graphical 
comparison is made of the two field descriptions for the parallel-plate waveguide. This is 
believed to be the first graphical display showing the interplay between the two alternative 
descriptions of the fields. Though the simple model of the parallel-plate waveguide was chosen 
because it best demonstrates the ray - mode equivalence, the results extend to other waveguide 
and cavity structures. It is demonstrated that for a waveguide separation large compared to 
wavelength, the ray and modal descriptions are in good agreement. 

Finally, the conclusions and some future work evolving from this thesis are discussed in 
Chapter 6. 
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CHAPTER 2 

MODAL ATTENUATION IN 
MULTILAYERED COATED WAVEGUIDES 


2.1 Introduction 

One method used in reducing the RCS of cavity structures is to coat the interior of the 
cavity with lossy dielectric or magnetic materials. Previous studies [13], [14] on the effects of 
coating involve modelling the cavity as a coated circular cylindrical waveguide terminated by a 
perfect electric conducting (PEC) short and studying the propagation and attenuation properties 
of the normal modes of the waveguide. There are also many other applications in which it is 
desirable to line the wall of a conventional waveguide with a layer of coating material [15]- 
[20]. The lining serves to alter the modal fields in the waveguide in order to achieve either 
more or less attenuation for certain modes. 

Pathak and Altintas [21] showed that the interior irradiation contribution to the RCS of a 
shorted waveguide is due mainly to a few low-order normal modes. By attenuating low-order 
modes, a lossy lining layer serves as a mode suppressor, which then reduces the RCS of the 
cavity. Lee et al. [14] studied the normal mode behavior for a circular waveguide lined with a 
single layer of coating. It was shown that at low frequencies (a/X ~ 1, where a = radius of the 
cylinder; X = the free-space wavelength) for a slightly lossy coating layer, the low-order 
modes can be highly attenuated with a single thin layer of coating. As the frequency increases, 
however, the attenuation constants of most of the low-order modes become small and decrease 
as a function of XV a?. Thus, RCS reduction can be achieved only over a fairly narrow 
frequency range with a single layer of coating. This chapter studies the propagation/attenuation 
constants of the low-order modes in a waveguide lined with multiple layers of coating. By 
using two or more layers of coating material, the fields inside the waveguide can be more easily 
manipulated to achieve a greater region of attenuation. With multilayered coating, it is possible 
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to achieve larger attenuation constants over a broader frequency range than with a single layer 
of coating. 

Past studies of coated waveguides have included such methods as perturbation theory 
[16], [17], [19], transmission-line model [15]-[18], and asymptotic theory [20]. These methods 
often require simplifying assumptions that are too restrictive, such as the coating material must 
be nearly lossless [16]-[18] or very lossy [20]. This treatment will apply the more general 
method of solving the modal characteristic equation exactly by a numerical method. This is 
feasible because of the fast computational speed of modem computers and the availability of 
efficient subroutines for computing Bessel functions with complex arguments. 

In Sec. 2.2, an overview of the modal fields in a coated circular waveguide is given. 
The mathematical formulation for the exact characteristic equation of the normal modes for a 
circular waveguide with multiple layers of internal coating is presented in Sec. 2.3. The 
conventional method involves setting the determinant of a 4n x 4n matrix equal to 0, where n is 
the number of coating layers. The propagation constants of the normal modes are the solutions 
of the characteristic equation. As the number of layers increases, this method becomes 
cumbersome due to the large size of the matrices involved. In addition, the characteristic 
equation, a transcendental equation involving Bessel functions, must be solved numerically on 
a computer. Different programs must be written for different numbers of coating layers. A 
generalized method for formulating the characteristic equation is presented which involves only 
the manipulation of 4 x 4 matrices. This method allows for a computer implementation which 
permits any arbitrary number of coating layers. Numerical results for the attenuation properties 
of the normal modes of circular waveguides with single and multiple layers of coating are 
presented in Sec. 2.4. 
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2.2 Overview of Modal Fields in a Coated Circular Waveguide 

A coated cylindrical waveguide is shown in Fig. 2.1 The problem of interest is the 
propagation/attenuation properties of the normal modes of such a structure. Figure 2.2 shows 
an exaggerated view of the coating layers for the sake of illustrating the geometrical features of 
the guide. In the application, the coating layers will be very thin relative to the diameter of the 
guide. The waveguide walls are assumed to be perfectly conducting, and each coating layer is 
assumed to be of uniform thickness. The axis of the cylinder coincides with the z axis. Both 
the permittivity ei and permeability |ii of each coating layer are allowed to be complex. The 
characteristic equation for the propagation constant kz of the normal modes is derived from the 
well-known method of seeking nontrivial solutions for the coefficients of the field expressions 
of the equations obtained by enforcing the continuity of the four tangential fields, Hz, Ez, H(|), 
E(j), at each interface between two coating layers, and between the innermost coating layer and 
the inner region [22]. 

In an uncoated waveguide, the normal modes are either TE^n or TMmn with respect to 
the longitudinal axis z. The index m describes the azimuthal variation in the form of sin m(j) or 
cos m<j). Here, the index n describes the order of the eigenvalues of Jm(kpa) = 0 for TM and 
J'm(kpa) = 0 for TE, where Jm is the Bessel function of order m. (Note that variable n is also 
used to represent the number of coating layers. The intended meaning of subsequent uses of n 
should be clear from its context.) With the exception of the m = 0 case, the normal modes are 
no longer pure TE or TM when the waveguide is coated with dielectric or magnetic material. 
The modes are commonly classified into HEmn and EHmn in such a way that in the limiting 
case of a vanishing thin coating [23] 


HEmn TEmn, 


and EHmn“^ TMmn- 


( 2 . 1 ) 
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The modal fields for the waveguide are given by 
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The convention of exp[j(cot - kz z)] is understood and suppressed. Superscript 1 
represents the inner region; superscript n+1 represents the layer of coating on the conducting 
wall; and superscript i represents the i^^ region, where 2 < i < n. Subscripts p, (|), and z 
indicate the radial, angular, and propagarion-directional components of the fields, respectively. 
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kp- represents the radial wave vector in region i, where k^pj + k\ = eijiik^Q and ko = 2n/X; (o 
is the angular frequency; Jm is the Bessel function and Nm is the Neumann function of order 
m. Ai, A 2 , BH, B^ 4 ^ Ci, and C 2 are the constants, which are determined by the 

boundary conditions and normalization requirements. Due to circular symmetry, there are two 
degenerate modes for each angular mode index except for m = 0. One of the two degenerate 
modes is arbitrarily chosen in the above expressions. 


2.3 Characteristic Equation of the Norma! Modes 

For the case of a single layer of coating, the characteristic equation for the propagation 
constant kz of a normal mode is well known [22] and given by 
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For a double layer of coating, the characteristic equation for the propagation constant kz 
of the modal fields is given by the 8 x 8 matrix equation 
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(* Ml 
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— — G,(b) 

(Ob Ml 
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m 1 
^ K 

(Ob M3 

k m 1 

— — Lf,« 

coa ^ 

k„ F.ta) 
Pi 1 

-km 1 
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(oa 

-k m 1 
' ' G,(a) 

(Oa 

-k„ F,ta) 
Pi 2 

-k G-ta) 
Pi 2 

0 

0 

0 

0 

LF,(b) 

(Ob £, ^ 

km 1 

—-Lgm 

(do £2 ^ 

K 

Pi 2 

Pi 2 

-km 1 

" * K,(b) 

(ob £3 ^ 

-kp^Kitb) 


= 0 


where a = ri, b = T 2 , and c =T 3 , (2.4) 

In general, for n layers of coating, the characteristic equation for kz is given by a 4n x 4n 
matrix equation. In numerically solving the equation, the computer time required becomes 
cumbersome as the number of layers increases. An efficient method for formulating the 
characteristic equation [24] for an arbitrary amount of layers can be utilized which involves 
only the manipulation of 4 x 4 matrices. 
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The boundary matching equations at each media interface can be written as 4 x 4 matrix 
equations. At p = ri (boundary between region 1, center region, and region 2, innermost 
layer). 


^Pi 

— Fi(ri) 


0 


kpFi'(ri) 


kzm 1 

OlTi £i 


0 


0 


•"Pi 

41 


Fi(ri) 0 


1 

COTi Pi 


Fi(ri) 0 


Fi(ri) kpFi'(ri) 0 


^2 


F2(ri) 


0 


k^m 1 


£2 


0 


kp/2'(ri) kpG/Cri) 


k,m 1 


0 


k*" 

^F,(n) ^G,(ri) 


k^m 1 k,m 1 

-F2(ri) -^7-G2(ri) 


(ori P2 


cori P2 


Off, 


-^i7G2(ri) kp F/Cri) kpG2'(ri) 


B 


B2* 




bJI 


(2.5) 


16 


At p = Tj, 2 < i < n - 1 (boundary between region i and i+ 1), 


fF,(r.) 

kp. 

^G.(,) 

0 

0 

B’l 

0 

0 

k^ 

^F,(„ 

*^pi 


kpFi'(ri) 

kpGi'(rj) 


' G.(,) 

WTi Hi ‘ ‘ 

B^ 

kzm 1 

k,m 1 

— -Gj(ri) 
corj £i * * 

kp_F/(ri) 

kp,Gi'(r,) 

b; 


ti+i 


0 


kztn 1 

L £i+i 


Pi 

Ei+l 


Gi+i(ri) 


0 


kpFi+i'(ri) kpGi+i'(rj) 


Pi' 


■Fi^i(ri) 


k5,m 1 

(OTi Ei+i 


0 




Gi+i(ri) kpFi+i'(r;) 


cori Pi+i 

kpGi+i'(ri) 






B 


i+l 


B 


i+l 


(2.6) 
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At p = r„ (boundary between two most outermost layers), 



I 



(2.7) 
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The boundary matching equations thus become n matrix equations for n layers of coating; 
Mil A = M21 B 
M22 B — M32 B 


Mnn B = M(„^1)„ C. ( 2 . 8 ) 

Mab refers to the matrix resulting from the tangential fields in region a matched at boundary rb. 
A is then related to C by 

Mil A = M21 M ^2 • • • M(n+i)„ C ( 2 . 9 ) 

or MiiA = MC 

where M = M21 . . . M^J, M(n+i)n. 

Because the last two columns of M(n+i)n are zero vectors, the last two columns of M will 
also be zero vectors. Therefore, M will be in form 


M = 


'mil 

mi2 

0 

O' 


m2i 

™22 

0 

0 

(2.10) 

msi 

^32 

0 

0 


mai 

m42 

0 

0 
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Equation (2.9) written out is then 


'‘p. 

0 

0 

0 



mil 

mi2 

0 

0 

Cl 

0 

^Pl 

TTfF,(r,) 

0 

0 

A 2 


m2i 

m22 

0 

0 

C 2 

kp,F'(r.) 

cor, 

0 

0 

0 


msi 

ms2 

0 

0 

0 

k,m 1 

L '"‘1 *^1 

k„.F/(r,) 

• 1 

0 

0 

0 


ni4i 

m42 

0 

0 

0 


which can be rewritten as 


r k^ 

-rFi(ri) 

0 

-mil 

-nii2 

^1 

0 


-m2i 

-m22 

^2 

kpF'(ri) 

k,m 1 

-mgi 

-m32 

Cl 

k,,m 1 

kpF,'(r,) 

-m4i 

-^42 

C 2 


(2.12) 


Ai, A 2 , and Ci, C 2 fully define the fields in the inner region and outermost layer, 
respectively. The characteristic equation for the propagation constant kz is then given by the 
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determinant of the 4 x 4 matrix Mac set equal to zero, where Mac is the 4 x 4 matrix in 

( 2 . 12 ). 


det Mac “ 


‘'p. 


k^m 1 


0 


Pi 

0 

H-l 


-mil “^12 


■^21 “™22 


k^m 1 

kpF'(ri) t;::-— -"^31 -^32 


COTi £i 


Fi(ri) kpFi'(ri) -m4i -m42 


= 0 


(2.13) 


To summarize, the boundary field matching equations can be expressed in the form of 4 x 
4 matrix equations (2.5),(2.6), and (2.7). From these relations, the coefficients for the fields 
in the center uncoated region can be shown to be related to the coefficients in the outermost 
coating layer by (2.9), which can be expressed as (2.12). Finally, for nontrivial solutions of 
the modal fields, the 4 x 4 matrix in (2.12) must be singular. Thus, the characteristic equation 
for the propagation constants of the normal modes is given by (2. 13). 

All manipulations now involve only the inversion and multiplication of 4 x 4 matrices. 
Using this formulation, a computer program to compute the propagation constant kz for any 
arbitrary number of coating layers can be written. The final expression (2.13), a complex 
transcendental equation, can be solved numerically, using, for example, Muller's method 
(available as an International Mathematical Statistical Libraries subroutine). 
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2.4 Numerical Results 

The attenuation constant a of a normal mode in a coated waveguide is related to the 
magnetic and electric energies in the coated regions, IHP and IE|2, by 

a « J^[Imlii|2|Hl2 + Imle|2|E|2] dV (2.14) 

where the volume of integration is over the coated region [22]. Because of the boundary 
condition of the PEC waveguide wall, the magnetic energy of an empty waveguide is much 
larger than the electric energy near the surface. Since the fields of the waveguide are not 
nerturbed sisnificantlv bv a coating laver in the low frennenrv recrion a loccy and magnetic 
coating material is effective in reducing the RCS of a coated waveguide. Lee and Lee [12] 
theoretically predicted and experimentally verified a 20 dBSM reduction for a PEC- terminated 
waveguide {afk = 0.98) lined with a lossy magnetic coating material only 1.18% radii thick. 
As such, the numerical results presented will focus on coating configurations involving 
magnetic layers. 

The attenuation properties of the two dominant low-order modes for near axial-incidence 
RCS, the HEii and EHn modes, will be discussed for single- and double-layer coating 
schemes. The interior uncoated region is taken to be free space (e = eo, |i = [to) all cases. 
The transverse field distributions for the two modes are shown in Fig. 2.3 for (a) empty guide, 
(b) lossless dielectric coated guide near cutoff, (c) lossless magnetic coated guide near cutoff, 
and (d) lossless coated guide in the high frequency limit. 

2.4.1 Review of single-layer coating results 

With a layer of lossy coating material, the low-order modes of a circular waveguide 
become inner modes [14] as aA increases. The field distribution for such modes is confined 
mostly in the center region. The fields decay very rapidly from the coating interface to the 
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Fig. 2.3 Transverse field distributions of the normal modes in a) empty waveguide b) dielectric 

coated waveguideguide (£r = 10.0, Pr = 1-0) at cutoff frequencies c) magnetic coated 

waveguide (£r = 1 .0, Pr = 10.0) at cutoff frequencies, and d) coated waveguide at the high 
frequency limit. 
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lossy layer, and the attenuation constants are small. Figures 2.4 - 2.6 show the attenuation 
constants as a function of increasing frequency for the HEn and EHn, as well as the HE 12 , 
modes in a lossy magnetic (jir =1.5 - j2.0) coated waveguide. As the frequency increases, the 
lossy material tends to expel the fields from the coating region. The attenuation constants 
decrease as a function of ^2/a3 for large frequencies [14]. 

If the coating material is lossless, however, the modal power distribution is largely 
concentrated in the coating region. A layer of high dielectric constant tends to "pull" the fields 
into the coating region as the frequency increases. As shown in Fig. 2.7, although the 
dielectric layer (er = 10.0) covers only 12 percent of the waveguide cross section, 99 percent of 
the power is confined in the dielectric layer. These results suggest that by using a double layer 
of coating consisting of a lossy layer and a lossless layer with a large permittivity, the high 
attenuation region of some of the modes of the waveguide can be extended. A high field 
concentration would be attracted by the lossless dielectric layer into the coating region, and 
attenuated by the lossy layer. In the following multilayered coating results, the lossless 
dielectric layer will have Er = 10.0, and the lossy magnetic layer will have pr = 1-5 - J2.0. 

2.4.2 Multilayered coating 

Figures 2.8 and 2.9 show the TEn attenuation constant as a function of increasing 
frequency for a lossless dielectric layer sandwiched between the waveguide wall and a lossy 
magnetic layer. Figure 2.8 shows a family of curves for a fixed thickness of lossy material and 
varying thicknesses of lossy dielectric material. As the dielectric layer increases in thickness, a 
greater amount of the field energy will be concentrated in this layer. This also causes more of 
the field to be within the lossy layer, resulting in greater attenuation. Note that there is a 
significant increase in both the level of attenuation and the frequency band of high attenuation. 
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Fig. 2.5 Attenuation constants of the EHi i mode in a lossy magnetic coated 
waveguide (£r= 1.0, |lr = 1.5 - j2.0). 
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Fig. 2.6 Attenuation constants of the HE 12 mode in a lossy magnetic coated 
waveguide (£r = 1.0, |J.r = 1-5 - j’2.0). 
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Fig. 2.7 Normalized angle-averaged power distribution in WattsA^ as a function 
of radial distance in a dielectric coated guide (£r = 10.0, Pr =1.0). 
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Fig. 2.8 Attenuation constants of the HEn mode in a double layer coated waveguide 
with an inner layer of lossy magnetic material (Er = 1.0, |J.r = 1-5 - j2.0) and 

an outer layer of lossless dielectric material (Er = 10.0, |ir = 1-0). The 
thickness of the inner layer is fixed while that for the outer layer is varied. 
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Fig. 2.9 Attenuation constants of the HE 11 mode in a double layer coated waveguide 

with an inner layer of lossy magnetic material (e, = 1.0, = 1.5 - j2.0) and 

an outer layer of lossless dielectric material (£r = 10.0, = 1.0). The 

thickness of the outer layer is fixed while that for the inner layer is varied. 
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In Fig. 2.9, the thickness of the lossy layer is varied while the thickness of the lossless layer is 
increased. If the lossy layer becomes too thick, the field repulsion property of the lossless 
layer is dominant over the field attraction property of the lossless layer. Figure 2. 10 reverses 
the lossless and the lossy layers so that the lossless layer is towards the waveguide center. No 
significant increase in attenuation is achieved with this configuration for the TEn mode. 

Finally, Fig. 2.11 shows the EHn attenuation constant as function of increasing 
frequency for the waveguide wall - lossless layer - lossy layer geometry. The family of curves 
represents a fixed thickness of the lossy layer and varying thicknesses of the lossless layer. 
Again, there is a significantly higher level of modal attenuation over a wider frequency band for 
the multilayered case compared to that of a single layer of the lossy coating case. Note that the 
maximum value of the attenuation axis is now 100 dB/a (a = radius of cylinder). It is 
interesting to observe that up to aA ~ 2.0, a thicker lossless layer results in a higher 
attenuation. But for aA > 2.0, the x = 0.01a case results in the highest attenuation. This 
geometry, though, is not an effective coating for the HEn mode (Fig. 2.8). 

2.4.3 Radar cross section reduction 

The application of interest is the reduction of RCS of coated cavity structures. In such 
applications, the actual coating geometry must be optimized to effectively attenuate all the 
highly excited modes. Detailed discussion of the evaluation of the RCS from a coated 
waveguide can be found in Reference [12]. Figure 2.12 shows the theoretical RCS 
calculations for a length of circular waveguide terminated with a flat PEC plate. The diameter 
of the guide is 4 X, and the length is 10 X. The three curves correspond to the RCS of (i) an 
empty guide, (ii) a waveguide coated with a single layer of lossy magnetic coating, and (iii) a 
multilayered coated waveguide consisting of a lossless dielectric layer sandwiched between the 
waveguide wall and a lossy layer. For the given size of the structure, the single layer coating is 
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Fig. 2. 10 Attenuation constants of the HEn mode in a double layer coated waveguide 
with an inner layer of lossless dielectric material (er = 10.0, |ir = 1.0) and an 

outer layer of lossy magnetic material (Er = 1.0, j^r = 1.5 - J2.0). The 
thickness of the inner layer is fixed while that for the outer layer is varied. 
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Fig. 2. 1 1 Attenuation constants of the EHn mode in a double layer coated waveguide 
with an inner layer of lossy magnetic material (8r = 1.0, jir = 1.5 - j2.0) and 

an outer layer of lossless dielectric material (£r = 10.0, |ir = 1-0). The 
thickness of the inner layer is fixed while that for the outer layer is varied. 
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not very effective in reducing the RCS. However, a double layer coating results in a 
substantial reduction in the RCS of the cavity at near axial incidence. When the angle of 
incidence of the incoming plane wave 0 = 20°, the three curves approach a similar level. This 
is due to the fact that the higher-order modes are now more strongly excited, and the chosen 
coating geometry is not effective in attenuating these modes. 

2.5 Summary 

This chapter studied interior coating as a means of reducing the RCS of arbitrary cavities 
by studying modal attenuation in coated waveguide structures. An overview of the modal 
fields in a coated circular waveguide was given. To avoid dealing with excessively large 
matrices, a method was shown to express the waveguide characteristic equation for any 
number of coating layers which involves only the manipulation of 4 x 4 matrices. This method 
was applied to study modal attenuation in a multilayered coated waveguide. 

The dominant contributors to the RCS of a PEC-terminated waveguide at near-axial 
incidence are the low-order modes. A thin single layer of lossy magnetic material can greatly 
attenuate these modes, and thus reduce the RCS of a waveguide cavity. However, this is only 
effective over a narrow frequency band at low frequencies. With the proper combination of 
lossless dielectric coating and lossy magnetic coating, it was shown that a significantly higher 
modal attenuation can be attained over a broader frequency band for the HEn and EHn 
modes. For large off-axis incidence, higher-order modes come into importance. It should be 
possible to reduce these modes over a broad frequency band with the proper combination of 
three or more layers of coating. 
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CHAPTER 3 

SHOOTING AND BOUNCING RAYS: 

CALCULATING THE RCS OF AN ARBITRARY CAVITY 

3.1 Introduction 

Traditionally, the problem of calculating the RCS of a cavity structure is treated by the 
modal analysis. For example, Chang and Senior [9] studied the open spherical shell problem 
based on the expansion of interior and exterior fields in the spherical wave function. Recently, 
the problem of scattering from an open-ended circular cylinder with wall coatings was analyzed 
by Lee et al. [25] by utilizing the cylindrical waveguide modes. Analytic solutions to a family 
of canonical problems have also been obtained by the dual series approach [26], [27]. 
However, the aforementioned techniques become impossible if (i) the shape of the cavity is not 
a perfect sphere or cylinder, (ii) the electrical dimension of the cavity is large, or (iii) the space 
inside the cavity is not homogeneous. These restrictions prevent the realistical modelling of 
physical problems and hinder the study of the effects of shaping on RCS. 

In this chapter, a different strategy for analyzing the open cavity problem, entitled 
shooting and bouncing rays (SBR), is presented. A dense grid of geometric optics (GO) rays 
representing an incident plane wave is "shot" into the cavity through the front aperture and 
followed as the rays bounce from conductors, penetrate through materials, and eventually 
return to the opening of the cavity (see Fig. 3.1). An innovative scheme is then used to 
integrate the aperture field to obtain the scattered field. This approach has the following 
features. First, this is the first attempt of its kind at solving partially open structures using the 
ray technique. This approach is different from the traditional ray-tracing or ray-mode hybrid 
methods. Second, a real physical problem can be modelled closely, taking into account the 
noncircular opening of the cavity, the wall coating, and the longitudinal bending or twisting of 
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Fig. 3.1 Ray bouncing approach to calculating the RCS of a partially open cavity 
with complex geometry and material loading. 
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the cavity. This method is so simple in concept that there is virtually no restriction on the 
shape or material loading of the cavity. 

In Sec. 3.2, the problem is formulated based on geometrical optics. The paths of each 
individual ray are first determined by Snell's law. The field amplitudes associated with each 
ray are computed by taking into consideration (i) geometrical divergence factor, (ii) 
polarization, and (iii) material loading of the cavity walls. The contributions to the 
backscattered field from individual rays are then summed up to arrive at the total RCS due to 
the interior irradiation of the cavity. Results of the SBR method are compared to the existing 
modal results for the case of a straight circular cylinder with a PEC termination in Sec. 3.4. 

3,2 Formulation of Shooting and Bouncing Rays 

Referring to Fig. 3.1, consider an arbitrary cavity with an opening at the aperture 

The inside walls of the cavity may be coated with dielectric or magnetic materials. Additional 
scatterers may also exist inside the cavity, but will not be considered in this treatment. The 
incident plane wave is given by (for exp(+jcot) time convention) 

Ii = [4l + §il]ej^^-r (3.1) 

where 

k i = ko (x sin 0^ cos ([>i + y sin 0* sin ())i + z cos 0^) 

I = amplitude of the TE wave 
I = amplitude of the TM wave. 
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The problem at hand is to determine the backscattered field in the direction of and the RCS 
of the cavity. The backscattering from the exterior of the cavity is not treated, since (i) the 
exterior scattering is small for small angle of incidence 0^ and (ii) the exterior scattered field can 
be considered separately. In addition, it is well known that the contribution to the 
backscattered field from the energy coupled into and re-radiated from the cavity is much 
greater than that from the rim diffraction, especially if the cavity opening is large compared to 
the wavelength [8], [10]-[12], [25], [28]. Therefore, attention will be focused on the interior 
irradiation mechanism. 

The ray bouncing approach will be carried out in three parts: 

(i) Given the geometry of the cavity and the incident field, find the ray paths in the cavity by 
ray tracing. This part of the problem is dependent only on the geometry of the problem. 

(ii) Determine the field amplitude of the exit rays on the aperture based on geometrical optics. 
This involves calculating the ray tube divergence factors and the reflection coefficients. 

(iii) Use Kirchhoffs approximation (physical optics) to determine the backscattered field and 
the RCS. 

These steps will be described in detail in the following sections. 

3.2.1 Ray tracing 

In order to model the incident plane wave, parallel rays are launched from the incident 
direction. Each ray is represented by a line in space with a reference point (xq, yQ, Zq) and a 
direction vector (sj, S2, S3). Any point (xj, y|, z^) along this line would then be described by 
(xj, y^, Zj) = (xq, yQ, Zq) -h (sp S2, S3) t . (3.2) 

By this definition, if the phase of the field at (xi, yi, zi) lags that of point (xq, yQ, Zq), then t 


will be a positive quantity. The direction vector of the incident rays is given by 
Sj = -sin 9^ cos (|)^ 


(3.3) 
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5 2 = -sin 01 sin ())i 

5 3 = -cos < 1)1 . 

The reference point (xq, Yq, Zq) on the incident plane '!} can be related to the point (xa, ya. 0) 
on the aperture Sa via 

xo = (S2 + S3) Xa - SjS2ya + s^tg 

yo = (sj -I- S3) ya - SiS 2 ya + s^tg (3.4) 

ZQ = -SiS3Xa - S2S3ya + S3tg , 


The parameter to in (3.4) determines how far is from the aperture plane and can be chosen 
arbitrarily (e.g., to = -10). To summarize, the incident ray is described by (3.2), where the 
direction vector is given in (3.3) and the reference point (xo, yo, zo) related to a 

corresponding point on the aperture via (3.4). 

Once the incident rays have been defined, the impact point of each ray on the inner wall 
of the cavity can be determined. This is accomplished by solving simultaneously the equation 
describing the cavity z = f(x,y) and (3.2). For example, if the cavity is a circular cylinder with 
radius a, the intersection is found by substituting into (3.2) 

t = (-B -t- Vb2-4AC) /2A (3.5) 

where 

A = s?-.s 2 

B = 2(SjXg + S 2 yo) 

C = Xg + yg - a2 . 
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For arbitrary cavity structures, the condition Ic i • n > 0, where n is the inward pointing normal 
of the wall, may be useful for discarding extraneous solutions. 

Next, the equation for the reflected ray will be determined. The reflected ray must satisfy 
Snell’s law, namely, (i) it must lie in the plane of incidence and (ii) the angle of reflection must 
equal the angle of incidence. Referring to Fig. 3.2, define a unit vector 

m = (]^ x n)/sin 0^ . (3.6) 


m is perpendicular to the plane of incidence. The above two conditions can be restated as 
(i) • m = 0 and (ii) • n = cos 0^. These two conditions are solved simultaneously in 

order to determine the direction of the reflected ray. The equation for the reflected ray is now 
complete since both the direction and the reference point (impact point) are known. 

Consequently, by using this reflected ray as the incident ray, the above procedure can be 
repeated until the ray exits the cavity. The exit of the ray from the cavity is easily detected 
when the ray intersects the aperture Xa- For every ray launched, a set of impact points inside 
the cavity and the direction of the existing ray are obtained. 

3.2.2 Amplitude tracking 

Once the ray paths inside the cavity are found, the field amplitude along the ray can be 
determined. For each ray launched, a set of impact points was found within the cavity, {(x^, 
y^, Z|)}, i = 1, 2, ..., N. The aperture field associated with the exit ray, ^(xj^, yj,^, 0), is to 
be determined. From geometrical optics, the electric field obeys the following recursion 
relationship 


^(xi+i,yi+i,zi+i) = (DF)i • (Hi • ^(xi,yi,zi) • e-j (Phase) 


(3.7) 
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where Phase = kg [(x^+j, - Xj)2 + (yj^j, - yj)2 + (zj+j, - Zj)2 ]l/2 and ^(xj ,y^ .z^) is the inci- 
dent field at (xj, y^, Zj). (F)i is the planar reflection coefficient matrix at the i^^ reflection point 
where the original curve interface in Fig. 3.2 is replaced by its tangent plane at the reflection 
point. (DF)i is the divergence factor which governs the spreading of the differential ray tube 
from just after the ith reflection to just before the (i+l)^h reflection. The remaining tasks are the 
determinations of these two quantities. 

A. Planar Reflection Coefficients. Referring to Fig. 3.3(a), attention will be restricted to 
the case of one layer of dielectric or magnetic coating with thickness x and backed by perfect 
conducting walls. The transmission line analogy to the field problem is shown in Fig. 3.3(b). 
The well-known planar reflection coefficients for the TE and the TM case are summarized in 
Table 3.1. 


Table 3. 1 Reflection Coefficients 


TE 

TM 

1 

1 

N 

0 

N 

1 

- 1 - Zo/Z-^ 

‘ ~ 1 + Zo/Z_> 

” 1 Zo/Z-^ 

Z-> = jZi tan (Pzix) 

Z-4 =jZi tan (pzix) 

Zi = C0|li/pzi 

Zi = pzi/coei 

Zo = O)|io/Pz0 

Zo = PzO/coeo 


pzo = cos 0* 


pzo = co^ 




2 0 


1 

c 
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n 




Fig. 3.3 (a) The original curved interface is replaced by its tangent plane at the 

impact point 1 for calculating the reflection coefficient based on GO. 
(b) Transmission line analogy of the planar reflection problem in (a). 



For the assumed time convention, the complex root is chosen to have a positve real 


component and a negative imaginary component. 

R. TE/TM Decomnosition . In order to apply the reflection coefficients F and F found in 
Table 3.1, the incident field needs to be decomposed into its TE and TM components. Using 
subscript c to denote local coordinates (see Fig. 3.2), the incident field can be written as 

I(l-) = (g‘ 4c)$c + (®‘-^c)^c ■ 


The reflected field is then given by 


t(l+) = T(P ■ ^l) + F(ti • . 


(3.9) 


So the reflected field is easily found once and 
coordinate. Choose 

A A A 

= m X n 

A A 

Zc = -n 

A A 

yc = -m . 

From Fig. 3.2, 


are determined in terms of the global xyz 


(3.10) 


9^ = cos'^ (^b • n) 

4 = 0 . 


(3.11) 
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Then 


^ cos 0^ - Zf. sin 6^ 


A1 A 

<t>c=yc 


and 


(3.12) 


cos 0^ - Zq sin 0^ 

Ar A 

<{>c = yc • (3.13) 


C. Divergence Factor . 
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3.2.3 Physical optics 

Given the outgoing field on the aperture the outgoing field can be computed by the 

standard physical-optics approximation. First, the outgoing field is replaced by an equivalent 
magnetic current sheet 


2 E(xn, y^, 0) X z ; over 


Ks 




I 0 


; outside 


(3.23) 


Kg radiates in the backscattering direction and gives rise to the RCS. The backscattered field is 
calculated from 


Ebs = l_[0i Aq + $i 


Ae 

f 

dxdye'’^'" 

cos <|)“ -f- Ey sin (])‘ 

_ A({)_ 

27t 

J 


(—Ex sin (|)‘ 4- Ey cos ({)*) cos 0* 


(3.24) 


u = sin 0^ cos (j)^ 

V = sin 0^ sin (])^ . 

and Ey are the x and y components of the outgoing field on the aperture E^. The RCS is 


defined in Table 3.2. 


preceding page 8LANK NOT FILMED 
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Table 3.2 Radar Cross Section Definition 



TE Incidence (I = 1, 1 = 0) 

TM Incidence (I = 0, 1 = 1) 

Co-pol 

RCS ^,4 = 4hIA/ 

RCS 00 = 4:cIA0|2 

Cross-pol 

RCS04=4itlA/ 

RCS^e = 4ltlA/ 

Phase 

^A^.ZAe 

^ A 0 ’ Z A^ 


In general, the positions of the outgoing rays are randomly dispersed over the aperture, 
even if the incident rays are launched uniformly. Since the outgoing rays do not lie on an 
equally spaced grid, the integration in (3.24) cannot be easily carried out. One possible 
solution to this problem would be to interpolate the results to find the fields on the integration 
grid. However, the proper interpolation scheme is difficult and an alternative approach is 
taken. Suppose only a small ray tube is shot into the cavity. This ray tube bounces around the 
walls and eventually comes to the aperture plane. It is possible to compute the backscattered 
field due to this ray tube by taking into account its wavefront curvature, size, and shape. By 
repeating this process until enough ray tubes are launched into the cavity to model the incident 
plane wave, the total scattered field should be the sum of the scattered field due to each 
individual ray tube. Using this idea, the scattered field expression will now be derived. 

Consider one of the incident ray tubes with an area of (Axq AyQ). The central ray with 
direction vector (sj, S 2 , S 3 ) hits point (xj, yj) on (see Fig. 3.4). The field within the 
exiting ray tube area (Axj Ay^) will be approximated as follows: 


"Ex(x,y)“ 


■E^(xi,yi)“ 

_Ey(x,y)_ 


_Ey(xi,yi)_ 


-jko[si(x - Xj) + S2(y - yi)] • 

6 


(3.25) 




Fig. 3.4 Planar approximation for the wavefront of the exit ray tube. 



In other words, the field with the ray tube at point (x, y) has the same magnitude as the field 
associated with the central ray. In addition, there is a linear phase variation across the ray tube. 
This linear phase approximation should be valid as long as the output ray tube is not too large. 
The physical-optics integral in (3.24) can now be evaluated using the plane-wave 
approximation. By summing over each ray tube, (3.24) becomes 


'V 

A, 




= ^ X J I dxdye 

2k i ith exit 


jko(ux+vy) -jko[si(x - Xj) + S 2 (y - yj)] 


exit 

all rays ray tube 


[ Ex(xi,yi) cos + Ey(xj,yj) sin ”| 

-Ejj(xj,y|) sin ())^ cos 0^ + Ey(xj,yj) cos (j)^ cos 0^J 


(3.26) 


Since the fields associated with each exit ray, Ej^(x|,yj) and Ey(xpyj), are independent of the 
integration variables x and y, the bracket in (3.26) can be taken out of the integral sign: 


■^0" 

" ^ S 

Ex(xi,yi) cos (|)^ + Ey(xj,yj) sin 


27t i 

O 7 1 f*0 x/c 

_-Ejj(xj,yj) sin (j)^ cos 0^ + Ey(xj,yj) cos cos 0^_ 


-jko[siXi + S2yil 


(AxjAyj) Ij 


(3.27) 
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h = 


1 

(AxjAyj) 


itn exit 


ray tube 


and 

(Ax^Ayj^) = area of the exit ray tube. 

Closer examination reveals that Ij is nothing but the Fourier transform of the ray tube 

shape function (normalized with respect to the ray tube area). In order to determine the shape 
of the ith exit ray tube, four adjacent rays are launched into the cavity in the form of a square as 

shown in Fig. 3.5. These four rays, upon exit, form a tetragonal shape on the ouqjui aperture. 
The position vectors of the four rays on the aperture are denoted by y = x^j x + yj^ y , n = 1 , 

2, 3, 4. The normalized Fourier transform of the tetragonal shaped function can now be 
evaluated as described in [34]: 

li = S(p,q) / S (0,0) (3.28) 

where 


S(va)= y ^ e“’ (^n+1 ~ ^nXYn ~ Vn-0 ~ (yn+1 ~ Xn)(^n ~ ^n-l) 

4 - Xn-l)P + (yn ' yn-l)qH(Xn+i - x„)p + (y^+j - yj^)q] 

S (0,0) = 1/2 I(xjy2 - X2yj) + (X2y3 - X3y2) + (X3y4 - X4y3) + (x4yj - Xiy4)l 


p = ko (u-Sj) 
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2 

Input Output 


Fig. 3.5 Shape of the exit ray tube on the output aperture is obtained by shooting 
four adjacent rays into the cavity. 
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q = kQ (V-S 2 ) 

-» A A 

w = px + qy . 

Note that by definition, Yq = 74 in the above expression. The shape of the ray tube is 
approximated by a four-sided polygon. If desired, higher-order polygons can always be used 
for a more refined approximation of the ray tube shape. 

Lastly, the output ray tube area (Ax^Ay^) can be related to the incident ray tube area 

(AxQAyQ) through the ray tube divergence factors calculated earlier. The area of the exit ray 
tube is related to the area of the incident ray tube via the product of the divergence factors: 

(AxiAy;) = [| l(DF)il] (AxoAy„) . ,3.29) 


The ratio of the cosines is an added factor which accounts for the angles of the ray tubes 
entering and exiting the cavity (see Fig. 3.6). 

If the quantity of interest is the bistatic cross section (BCS) rather than the assumed 
monostatic RCS, the incident direction angles 0* and ()>* should be replaced by the observation 
direction angles 0°*’ and in (3.24), (3.26), and (3.27). 

To summarize, the physical-optics integral defined in (3.24) is evaluated by summing the 
backscattering contributions due to each ray tube. The wavefront of the exiting ray tube is 
approximated to be planar as given in (3.25). The shape of the exit ray tube is approximated 
by a four-sided polygon. The far field is proportional to the Fourier transfomi of this 
normalized shape function and is found using (3.28). Finally, the area of the exit ray tube is 
found by the product of the divergence factors and is related to the incident ray tube via (3.29). 
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3.3 Numerical Results 

A ray bouncing program, denoted as SBR (^hooting and Bouncing Rays), has been 
implemented based on the formulation described in Sec. 3.2. In this program, it is assumed 
that the cavity is composed of an arbitrary cylindrical cavity with a flat aperture at z=0 and a flat 
PEC termination at z=d. The iimer walls of the cylinder can be coated with either dielectric or 
magnetic materials. Below, the comparison of SBR results with those generated by the modal 
analysis for a circular cylinder with uniform cross section is shown. 

3.3.1 Exit ray positions 

Shown in Fig. 3.7 are plots of the incident and exit ray positions on the aperture plane for 
a circular cylinder with a flat end-plate. The incident rays are launched from a uniformly 
spaced grid on the aperture. For a small incident angle (0 = T), the majority of the incident 
rays simply reflect off the end-plate and bounce back to the aperture without hitting the side 
walls. Therefore, they remain uniformly spaced on the exit aperture. As the incident angle 
increases, the positions of the exit rays become more dispersed due to the large number of 
bounces inside the cavity. This behavior is even more apparent in the case of an elliptical 
cylinder with nonuniform cross sections as shown in Fig. 3.8. 

3.3.2 Comparison with modal analysis 

RCS results generated by the SBR program are compared with those obtained by the 
modal analysis [12] for a straight circular cylinder with PEC walls. Shown in Figs. 3.9-3. 1 1 
are RCS plots versus angle 0 for circular cylinders with a fixed diameter of 4 X and a different 
depth d. The copolarization of both the TE and TM cases are plotted as RCS(j)(j) and RCS0 0 , 
respectively. The solid circles are the modal results and the solid curves represent results from 
the present SBR formulation. Good agreement is obtained between the ray bouncing method 




Fig. 3.7 The incident and exit ray positions on the aperture for a circular cylinder 
with a flat end-plate. 











RCS00(dBX^) RCS^^^CdBX ) 


(a) (])(})-polarization 
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Fig. 3.9 Comparison of results for ray bouncing with modal analysis for a 

circular cylinder with diameter 4X and depth 2X. (a) <})(l)-polarization, 
(b) 00 -polarization. 



(a) <l)(t)-polarization 
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Q (DEGREES) 

Fig. 3. 10 Comparison of results for ray bouncing with modal analysis for a 

circular cylinder with diameter 4^ and depth AX. (a) (|)(J)-polarization, 
(b) 0 0 -polarization. 



RCSee (dB X^) RCS,^^ (dB X^) 


(a) ([xt)-polarization 
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9 (DEGREES) 

Fig. 3.11 Comparison of results for ray bouncing with modal analysis for a 

circular cylinder with diameter 4X and depth lOX,. (a) <txl)-polarization, 
(b) 9 0 -polarization. 
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and the modal analysis, except in Fig. 3.1 1(a) for the 10 A, deep cylinder. Possible reasons for 
the discrepancy will be addressed in Chapter 5. Although a comparison with a much larger 
cylinder is preferable, the modal analysis program was limited to about a 4 X, opening. 

The convergence of the present ray formulation is checked by increasing the total number 
of rays launched into the cavity. Shown in Fig. 3.12(a) is the RCSee at 0 = 25” plotted as a 
function of the number of rays per wavelength. It is observed that the results converged for a 
density greater than about 10-15 rays per wavelength. As the depth of the cylinder increases, 
the rate of convergence slows down. This is presumably due to an increased number of ray 
bounces inside the cavity and the highly divergent exit ray directions on the output aperture. In 
Fig. 3.12(b), the RCS 99 at 0 = 50' is plotted as a function of the number of rays. At this 

larger angle, convergence is slower compared to the case in Fig. 3.12(a), again due to the 
increased number of bounces. 

3.3.3 Radar cross section of large circular cylinders 

Because the number of propagating modes increases exponentially with the radius of the 
cylinder, large aperture problems are not easily handled by the m.odal analysis. Modal analysis 
quickly becomes extremely cumbersome for large apertures because the propagation constant 
of each normal mode must be known. Since GO is a high frequency approximation, SBR is 
ideally suited for large aperture problems. Shown in Figs. 3.13 and 3.14 are the RCS plots of 
large circular cylinders. In Fig. 3.13, the cylinder has a diameter of 5X. and depth of lOX,. In 
Fig. 3. 14, the cylinder has a diameter of lOX and depth of 40X. Fifteen rays per wavelength as 
suggested by the convergence study were used. The results were generated on the NSF Cray 
X-MP supercomputer at the University of Illinois. 


RCSee (dBX^) RCS^<^(dBX^) 


Fig. 3.12 


(a) 64 




NO. OF RAYS PER WAVELENGTH 

Convergence of the RCS results as the ray density launch is increased, 
(a) e = 25°, (b) e = 50°. 
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3.4 Summary 

A ray bouncing method was presented for calculating the electromagnetic scattering from 
an arbitrary open cavity. This approach is based on tracking a large number of rays launched 
into the cavity through the opening and determining the geometrical optics field associated with 
each ray by taking into consideration the (i) geometrical divergence factor, (ii) polarization, and 
(iii) material loading of the cavity walls. The contributions to the backscattered field from 
individual rays launched into the cavity are summed up to obtain the total internal RCS of the 
cavity. It was shown that the RCS results obtained by the ray approach agree well with those 
for the modal analysis even for a small 4 X cavity. 
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CHAPTER 4 

RADAR CROSS SECTION REDUCTION 
OF CAVITIES BY SHAPING 


4.1 Introduction 

In addition to coating with lossy or absorbing materials, another method for reducing the 
RCS of a target is the proper shaping of the target . Through the use of shaping and radar 
absorbing materials, the RCS of a B-1 bomber is two orders of magnitude less than that of a 
B-52 [35]. A target should not have any flat, cylindrical, or conical surfaces which might be 
illuminated by the incident radar field along the normal to the target. An example of a target 
shape that yields a low RCS value over a wide range of incidence angles is the cone-sphere [1]. 
This chapter examines the effects of shaping, specifically, the effects of longitudinal bending, 
on the RCS of cylindrical cavities by applying the SBR method developed in Chapter 3. 

For a large (relative to wavelength) flat-plate terminated cylinder, there is a large RCS 
value at near axial incidence due to the reflection from the termination (for examples, see Figs. 
3d 3 and 3.14). The near axial RCS of such structures can be reduced by a smooth gradual 
bending of the cylindrical axis. One such longitudinal bend, an S-bend, can be described by a 
half cosine function. For a cylindrical cavity with the cavity aperture on the x-y plane, the 
body restricted to the z < 0 half-space, and the cylindrical axis coincident with the z-axis (e.g., 
Fig. 3.1), the axis of the cylinder is deformed to 

(xc, yc) = (0., H/2 [cos(7Cz/d)] - 1] ) for 0 < z < -d (4.1) 

where H is the amount of offset and d is the depth of the cavity. If H is greater than the largest 
radii of the aperture, no GO ray at axial incidence will directly reflect from the termination 
plate. 
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Two cavity geometries and their corresponding S-bend geometries will be considered. 
First is a circular cylinder terminated by a flat shorting plate. The second is a cylindrical 
waveguide with a triangular aperture and a circular end-plate. The transition extending from 
the triangular cross section to the circular cross section is smooth. The SBR method is 
currently the only available viable method for computing the RCS for such geometries. Section 
4.2 will describe the mathematics of the cavities and their SBR implementation. Numerical 
results will be presented in Sec. 4.3. 

4.2 Cavity Geometries 

In this section, the mathematical descriptions of the cavities to be analyzed will be 
presented. Section 4.4.2 describes the straight and S-bend circular cylinders. The formulation 
for the triangular-to-circular waveguide is given in Sec. 4.4.3. 

4.2.1 Circular cylinder 

The terminated circular cylinder and S-bend circular cylinder are shown in Fig. 4.1. The 
surface for thestraight cylinder is described by 

x2 + y 2 = for 0 < z < -d (4.2) 

x2 + y 2 < a2 for z =- d 

where a is the radius and d is the depth of the cylinder . The equation for a terminated S-bend 
circular cylinder is given by 

for 0 < z < -d 


x2 -H (y - H/2 [cos(7Cz/d)] - 1])2 = a2 
x2 -I- (y - H )2 < a2 


for z =- d 


(4.3) 
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where H is the amount of offset of the end-plate with respect to the aperture. If H = a, the 
offset of the end-plate is exacdy one diameter. When H = 0, (4.3) reduces to (4.2). Note that 
for the S-bend cylinder, a cross-sectional view at any z = constant plane also yields a perfect 
circle. 

An impact point (xi, yi, zi) of a ray with the cavity is found from solving the equation of 
the cavity parametrically with the equation of the ray 

(x, y, z) = (XQ, yo, zq) + (si, S 2 , S 3 ) t (4.4) 

where 

(xQ, yo, zq) is the initial starting point or the previous reflection point or the ray, 

(si, S 2 , S 3 ) is the direction vector of the ray. 

For the straight circular cylinder, the solution can be given in closed form by 

t = (-B + Vb2 - 4AC) /2A (4.5) 

where 

9 9 

A = S] + S 2 
B = 2(siXq + S2 yo) 

C = xJ + yJ-a 2 . 

The coordinate (xi, yi, zi) is then found by substituting (4.5) into (4.4). For the S-bend 
circular cylinder, the simultaneous solution to (4.3) and (4.4) is found numerically. 

In the computation of the curvature matrix, the obvious choice for the (u, v) coordinates 
of the straight cylinder is (u, v) = (<|), z), where (j) is the standard angular coordinate in a 
cylindrical coordinate system. For the end-plate at z=-d, (u, v) is set to (u, v) = (x,y). The 
cavity can then be described by 


for 0 < z < -d 
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r = (x, y, z) = (a cos (j), a sin ({), z) 
r = (x, y, -d) with x^ + y2 < a 


for z = -d 


The curvature matrix for the cavity, defined by (3.18), can be easily shown to be 


r 1/R 0 -I 

L 0 0 J 


for 0 ^ z < -d 


r 0 0 1 

L 0 0 J 


for z = -d 


For the S-bend, define an angle (p such that 


X = a cos(p 

y - H/2 [cos(jtz/d)] - 1] = a sin tp 


In other words, cp can be described by 


(p = tan* 1 {(y- H/2[cos(7Cz/d)]- l])/x} 


In this case, the choice for (u,v) is now ((p, z), and 


r = (x, y, z) = (a cos (p, a sin (p + H/2 [cos(7iz/d)], z) for 0 < z < -d 


(4.10) 


r = (x, y, -d) with x^ + y2 < a 


for z = -d 


The curvature matrix can be found from (3.18). 


4.2.2 Triangle-to-cirde transition waveguide 

A segment of the line on the z = 0 plane described by y = mx + b can be smoothly 
transformed in the -z direction to an arc of the circle described by x^ + y2 = a^ on the z= -d 
plane with the following equation: 
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a z (x2 + y2 - a2) + (z + d) (y - mx -b) = 0 for 0 < z < -d (4. 1 1) 


where a is a transition parameter. This can also be represented in cylindrical coordinates by 

p(((),z) = -g(z)(sin<l)-msin(J)) + 'y[g(z)(sin (>- m sin (|))]2+ a^+ 2g(z)b (4.12) 
z + d 

where g(z) for 0 < z < -d 

2az 


Figure 4.2 shows such a transition for a = 0.05, d = 45, m = 0, b = -16/V3, and a = 15, 
where ^ is restricted to -150° < <() < -3(P. Several of these surfaces can be combined to form a 
smooth transition from an arbitrary polygon to a circle. Care must be taken to insure continuity 
of the structure at the surface junctions. This is accomplished with the proper choice of a . 

Figure 4.3 shows such a waveguide structure with a triangle-to-circle transition. For a 
triangular aperture with side s, center at (xc, yc) = (0., 0.), and base parallel to the x-axis, the 
parameters that define the waveguide are 

ai = ao mi = 0.0 bi = -s/(2V3) -150o < (]) < -30° 

a 2 = 2ao m 2 = -V3 bi= s/-/3 -30o<<|)<90o (4.13) 

tt3 = 2ao m 3 = V3 b 2 = s/V3 90° <<))< -150o • 


In Fig. 4.4, cross-sectional cuts at various lengths along the waveguide for «o = 0.05 are 
shown. 

The equation for the S-bend transition is 

a z (x2 + (y - H/2 [cos(7tz/d)])2 - a2) + (z + d) (y - H/2 [cos(7cz/d)] - mx -b) = 0 


for 0 < z < -d . 


(4.14) 




Fig. 4.2 Surface with transition from line segment to arc of circle. Circle is 
centered about (x,y) = (0,0). 
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Fig. 4.3 Triangle-to-circle transition vv'aveguide. 
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Fig. 4.4 Cross-sectional view of transition waveguide at various distances along 
the axial direction. 
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Figure 4.5 shows the surface shown in Fig. 4.2 with the addition of an offset H = 27. Figure 
4.6 shows the S-bend triangle-to-circle transition waveguide. In computing the surface 
curvature matrix, define (p from 

X = p((p, z)costp (4.15) 

y - H/2 [cos(7tz/d) - 1] = p((p, z) sin (p 

where p((p, z) is given by (4.12). Thus, (p is exactly equivalent to (4.9). The choice of (u,v) 
coordinates, as for the S-bend circular cylinder, is (u, v) = (rp, z). Ths surface given 
parametrically by (rp, z) is then 

— > 

r = (x, y, z) = ( p(cp, z)cos tp, p(tp, z) sin (p + H/2 [cos(:cz/d)], z) 

forO:Sz<-d . (4.16) 

Equations (4.14), (4.15), and (4.16) reduce to a straight transition waveguide when H = 0. 
The intersection points for the rays with the cavity is determined from the solution of 
simultaneously solving (4.4) and (4.14). The equations for the flat end-plate at d = 0 are 
identical to the second half of (4.2), (4.3), (4.6), and (4.7) in Sec. 4.2.2. 

4.3 Numerical Results 

RCS computation results will be presented for the end-plate terminated straight and S-bend 
circular waveguides and the straight and S-bend triangle-to-circle transition waveguides. For 
the circular waveguides, the diameter used is 5X across and the length along the z -axis is 10^. 
The diameter offset parameters H for the circular S-bend are 2.5X and 5X,, so that the back 
plate of the cylinder is exactly one half and one diameter below the aperture, respectively. For 
the transition waveguide, the length of one side of the triangular aperture is SAX and the 
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diameter of the end circle is 4.5X. The length along the z-axis is 6.75^, and the offset for the 
S-bend case is 4.05^. A ray entering the cavity at axial incidence will not directly hit the 
circular end-plate. 

4.3.1 Exit points and ray paths 

Shown in Fig. 4.7 are the exit ray pictures on the aperture plane for the circular straight 
and S-bend cavities. Each dot represents the point on the aperture plane where a ray exits the 
cavity. The incident rays are launched from a uniformly spaced grid on the aperture. At axial 
incidence (0 = (j) = 0°), the incident rays reflect straight back from the end plate, and the exit 
ray positions are coincident with the input ray positions on the aperture plane. For the S-bend 
geometry, because there is no direct reflection, the output rays are dispersed even at axial 
incidence. Figure 4.8 shows some exit ray pictures for the straight and S-bend transition 
waveguides. 

Figures 4.9 and 4.10 show the ray paths for a single individual ray for the S-bend circular 
cavity. The figures show the ray path from three different perspectives. Each dot represents a 
reflection point inside the cavity. The area ratio given on each figure is the ratio of output area 
to input area for the differential ray tube associated with the ray, calculated from (3.29). A 
typical ray can take anywhere from a few bounces to well over a hundred bounces for the given 
geometry. The ray in Fig. 4.9 takes 5 bounces inside the cavity while the ray in Fig. 4.10 
takes 123 bounces inside the cavity. Both rays represent a plane wave incoming at axial 
incidence. The ray in Fig. 4.10 is similar to a ray for the whisper gallery mode [36]. SBR is 
ideally suited for computer-aided design (CAD) systems in the sense that a high RCS in a 
particular aspect angle can be traced back to the contribution from a particular portion of the 
target through examination of the rays paths and field amplitudes. The particular portion of the 
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Fig. 4.8 Exit ray pictures for (a) straight triangle-to-circle transition waveguide with flat 
plate end-plate and (b) S-bend straight triangle-to-circle transition waveguide 
with flat end-plate. 
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Fig. 4. 10 Ray path in an S-bend circular cylinder (5X offset). 


85 


target can then be reshaped or coated with absorbing material, and the RCS of the modified 
geometry can be recomputed. 

Figure 4.11 shows the front aperture for the S-bend triangle-to-circle transition 
waveguide. Each dot represents an input position for a ray to be launched at axial incidence. 
Figures 4.12 - 4.17 show the individual ray paths. The number next to each dot in Fig. 4.11 
indicates the number of the figure that shows the ray path for that particular input point. Again, 
the area ratio for each individual ray is given. Note for this sampling of six rays, the value of 
the area ratio varies from a low of 3.52 (Fig. 4.13) to a high of 477558 (Fig. 4.15). For 
highly curved closed cavities, a ray that takes a very high number of bounces generally has a 
very high area ratio due to a very low divergence factor product. No conclusion can be reached 
for rays that take only a few bounces. The ray in Fig. 4.12 has a higher area ratio than the ray 
in Fig. 4.16 despite making more bounces. 

4.3.2 Radar cross section computations 

Convergence of SBR for the S-bend circular cylinder is checked by increasing the number 
of rays launched into the cylinder. Figure 4.18 shows the BCS calculations for the 2.5 X 
offset S-bend cavity for an incident plane wave at (0i, <t>i ) = (10°, -90°). Figure 4.19 shows 
the BCS calculations for the 5.0 X offset case. A greater cavity bend requires a higher density 
of rays for convergent BCS calculations. Figures 4.20 and 4.21 show the BCS's for the 
straight and S-bend transition waveguides, respectively, for the same incident plane wave. 
Note that the S-bend transition waveguide requires twice the linear density for convergence 
over the straight waveguide. This results in a four-fold increase in computation time due to a 
corresponding increase in the number of rays. Figure 4.22 compares the calculated RCS for a 
2.5^ offset S-bend cylinder against a straight circular cylinder. Figure 4.23 compares the 
calculated RCS for a 5.0X offset S-bend cylinder against a straight circular cylinder. Though 
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Fig. 4. 1 1 Triangular aperture of transition waveguide. Each point corresponds to input 
point of ray for Fig. 4.12 - Fig. 4.17. 




89 



Fig. 4. 14 Ray path in an S-bend triangle-to-circle transition waveguide. 
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Fig. 4.15 Ray path in an S-bend triangle-to-circle transition waveguide. 
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(a) <|)(l)-polarization 



93 


(b) 00 -polarization 



OBSERVATION Q (deg) 


Fig. 4. 1 8 BCS with different input ray densities for an S-bend cylinder with diameter of 
5X, depth of lOA, , and offset of 1.5X. Direction of incidence 
(0in, <t>in) = (10° -90°). (a) <))<l)-polarization, (b) 00-polarization. 



(a) (t)(t>-polarization 



(b) 00-polarization 



Fig. 4. 19 DCS with different input ray densities for an S-bend cylinder with diameter of 
5^, depth of lOX, and offset of 5.0^. Direction of incidence 
(0in> 0in) = (10°, -90°). (a) <))<t)-polarization, (b) 00 -polarization. 




(a) (t)({)-polarization 



(b) 00 -polarization 



Fig. 4.20 BCS with different input ray densities for a triangle-to-circle transition 

waveguide with triangle side of 5.4X, end circle diameter of Fig. 4.5^, and 
depth of 6.75A,. Direction of incidence (0in, <{>in) = (10o -900). 

(a) ({)(l)-polarization, (b) 00-polarization. 
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(a) ([)(|>-polarization 



(b) 09 -polarization 



Fig. 4.21 BCS with different input ray densities for an S-bend triangle-to-circle transition 
waveguide with triangle side of 5.4>«, end circle diameter of 4.5X, depth of 
6.75X, and offset of 4.05X,. Direction of incidence (Gin, <t>in) = (10°, -90°). 

(a) <|)<|)-polarization, (b) 00-polarization . 





RCS^a(dBX ) RCS, .(dBX'^) 





(b) 00 -polarization 



0(deg) 


Fig. 4.22 Comparison of RCS's for straight circular cylinder (diameter of 5X and depth of 
lOA.) with 2.5 A, offset S-bend cylinder (same diameter and depth). 

(a) (l)<))-polarization, (b) 00 -polarization. 




Fig. 4.23 


Comparison of RCS's for straight circular cylinder (diameter of 5A and depth of 
10^) with 5.0X, offset S-bend cylinder (same diameter and depth). 

(a) dxb-polarization, (b) 6 9 -polarization 
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both S-bend geometries show a reduction of RCS, there is a significant higher reduction near 
axial incidence by the 5.0X, offset. As the incidence angle increases to ±20°, the RCS levels of 
both offsets become comparable. Figure 4.24 compares the RCS of a coated S-bend cylinder 
(2.5X, offset, 6r = 2.5 - jl.25, = 1-6 - jO.8) to that of the uncoated S-bend and the straight 

cylinders. The coated S-bend case shows a 30 dB - 40 dB reduction across the range of angles 
calculated. Finally, Fig. 4.25 compares the RCS of the straight, the S-bend, and the S-bend 
with coating (Er = 2.5 - jl.25, jir = 1-6 - jO.8) transition waveguides. Again, the uncoated S- 
bend shows a good RCS reduction at near axial incidence, but the coated S-bend shows an 
even greater reduction over a greater range of incident angles. 

4.4 Summary 

The bouncing ray method presented in Chapter 3 was applied to study the electromagnetic 
scattering from several cavity structures. The emphasis was on shaping as a method of 
reducing the RCS of the cavities. It was shown that with a longitudinal S-bend, the RCS's of 
cylindrical cavities can be significantly reduced at near axial incidence. However, an S-bend 
cavity with coating causes an even greater RCS reduction, and over a greater range of incident 
angles. Convergence of the SBR method was demonstrated. S-bend geometries require a 
higher density of rays for convergence than their straight cylindrical counterparts. 

SBR is a computer intensive method. The results in this chapter were generated on the 
NSF Cray X-MP supercomputer at the University of Illinois. For an S-bend structure, 
approximately 90 percent of the computation time is spent finding the ray paths in the cavity. 
In applications of this method, it is often of interest to compute the RCS of a geometry over a 
band of frequencies. The paths of the rays and the divergence factor (3.21) are independent of 
frequency. If the RCS is desired for more than one frequency value, the impact points of the 


100 


(a) 4xt>-polarization 



(b) 00-polarization 



Fig. 4.24 Comparison of RCS's for straight circular cylinder (diameter of 5X, and depth of 
lOX.) with 2.5X, offset S-bend cylinder (same diameter and depth) and same 
S-bend with coating layer ( thickness = 0.1^, £r = 2.5 - jl.25, Pr = 1-6 - jO.8). 
(a) <jxJ)-polarization, (b) 00 -polarization. 
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(a) <l)<t)-polarization 



(b) 00 -polarization 



Fig. 4.25 Comparison of RCS's for triangle-to-circle transition waveguide (triangle side 

of 5^, circle diamter of 4.5, and depth of 4.05X) with 4.05^ offset S-bend 
transition waveguide (same aperture, end diameter, and depth) and same 

S-bend with coating layer (thickness = 0.1 X, £r= 2.5 - jl.25, Pr= 1-6 - jO.8). 
(a) <t)(l)-polarization, (b) 00 -polarization 
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rays and the corresponding divergence factors could be written to a file and saved on the 
computer and simply read back for different frequency values. The density of rays used 
should correspond to a sufficiently high number to insure convergence for the highest 
frequency. This would result in a substantial reduction in total computation time. 
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CHAPTER 5 

RAYS VERSUS MODES: PICTORIAL DISPLAY OF ENERGY FLOW 
IN AN OPEN-ENDED WAVEGUIDE 


5.1 Introduction 

Although the SBR method has shown promise in treating complicated structures, one of 
the important issues to address is, What are the limitations of this approach? This chapter 
addresses the above question by considering the problem of a plane wave impinging on a semi- 
infinite parallel-plate waveguide. Two alternative descriptions for the fields propagating in the 
waveguide are the "exact" modal analysis and GO - ray analysis, used by the SBR method. A 
ccmpanscn is made of the two field descripLluii^ for the paraiiel-piate waveguide. Though the 


simple model of the parallel-plate waveguide was chosen because it best demonstrates the ray - 
mode equivalence, the results extend to other waveguide and cavity structures. 

Based on ray optics, the fields and, therefore, the energy flow inside the guide should 
exhibit a beam behavior shown in Fig. 5.1. This physically intuitive ray picture is shown to be 
valid to a good extent if the separation of the waveguide plates is large compared to the 
wavelength. However, as the waveguide opening becomes smaller, the beam picture begins to 
blur after propagating a certain distance into the guide. These results are consistent with the 
fact that GO fields are high-frequency asymsptotic solutions to Maxwell’s equations. Although 
the complementary roles of ray and mode have been studied extensively in the past [37] -[39], 
this is believed to be the first graphical display showing the interplay between the two 
alternative descriptions of the fields. 


5.2 Formulation 

The problem under consideration is sketched in Fig. 5.2. The waveguide separation is 
given by a and the angle of the incident plane wave is 9. For the transverse electric (TE) 
polarized case, the incident electric field is assumed to be (exp (jeot) time convention) 
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^^Ag-Jko(sin0x + cos0z) 

where Icq = ayyfiie is the free-space wavenumber. The field excited inside of the waveguide 
can be described by TE modes: 


„ “ „ . njcx -jk^z 

Ey - S sin -^e 
n=l 


Hx = I, — 

n=l \oyii 


[ . nTCx -jk^z 
sm— e-' ^ 


oo 



. nTC \ nTCx -jk_z 

1 cos— — -e ^ 

J ^ 


(5.2) 


Here k^ = [kQ - (mt/a)^] and Cjj is the modal coefficient of the «th mode. The sign of k^ is 

to be taken as either positive real or negative imaginary. 

For waveguide separation kQa » 1, Kirchoffs approximation can be used to determine 

the modal excitation coefficients by matching the incident field in (5.1) to the modal field in 
(5.2) at the aperture z = 0. The results are 
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^ 2 , ^ -jkn sin 0x . nTCx , 

Cn =7 I e ^ sin— dx 


(mr/a) 




- (^0 




sin G)" 


nTC 




if Icq sin 0 ~ • 

. ^ « njr 

if Kq sin 0 = ~ 


(5.3) 


It is important to point out that the present interest is mainly in observing the field behavior 
deep into the waveguide, not the details of the field distribution near the wave opening. 
Therefore, although the exact Wiener-Hopf solution exists for this problem [40], the 
Kirchhoff s approximation is sufficient for our purpose. 

The above described procedure can also be applied to the transverse magnetic (TM) case. 
For the sake of brevity, the results for the TM case are not included here. The conclusions 
reached for the TE case can also be drawn for the TM case. 


5.3 Numerical Discussion 

Shown in Figs. 5.3(a) and 5.3(b) are plots of the energy flow inside the guide for 
waveguide separation a = 50 X. and plane wave incident at 0 = 30°. Each arrow contains 
information on both the magnitude and the direction of the Poynting vector Re (^ x fl*) at that 
point. The size of the arrow is proportional to the magnitude of the Poynting vector. The 
direction of the arrow is the direction of the energy flow. All vectors with a magnitude of less 
than 3 dB of the incident field are too small to appear. The physically intuitive picture based on 
the ray argument is shown in Fig. 5.3(a) and the actual energy flow computed from the modal 
sum in (5.2) is shown in Fig. 5.3(b). Figure 5.3(c) is an enlargement of a section of Fig. 


RAY PICTURE 



Fig. 5.3 Plots of Poynting vectors inside of a 50A, waveguide based on (a) the 
intuitive ray picture, (b) the actual fields obtained by summing up 
modes, and (c) an enlargement of a section of (b). 
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5.3(b) to show the details of the Poynting vectors. The evanescent modes are of little 
importance in the modal sum in the sense that there is no noticeable change in the energy plot if 
they are included or discarded. The resemblance of the actual distribution to the ray picture is 
apparent in this case. The fields in Fig. 5.3(b) are the sums of 99 propagating modes with 
different phase velocities. It is fascinating to observe that the modes add up to give an orderly 
ray picture. 

As the waveguide separation is decreased to 10 and 3 X shown respectively in Figs. 5.4 
and 5.5, the actual field no longer retains the beam behavior after some distance into the guide. 
The blurring occurs at approximately 70 X into the guide in Fig. 5.4(b) and 5 X, in Fig. 5.5(b). 
In other words, the simple ray-optical picture is not accurate for describing the field after some 
depth into a small waveguide. Shown in Fig. 5.6 are the energy flow plots for various incident 
angles of a 5 X waveguide. The beam blurring is not a very strong function of the incident 
angle. 

Finally, an interesting phenomenon is observed in Fig. 5.5(b) where the beam behavior 
reemerges after many wavelengths into the guide. At approximately z = 25.5 X, the condition 

[Icq - (ntc/a )] = integer multiple of 2k (5.4) 


is almost satisfied for all five of the propagating modes. The result is that the pattern from 0 to 
25.5 X will be repeated approximately. Note, however, that the reemerged beam pattern no 
longer corresponds to the ray picture in Fig. 5.5(a). The same phenomenon is displayed in 
Fig. 5.7 for a 5 X waveguide. In this case the condition in (5.4) is satisfied at around z = 60 X 
for the nine propagating modes. If the waveguide separation is large and many propagating 
modes are excited, the distance inside the waveguide before observing this repetitive behavior 
will be extremely long. 
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intuitive ray picture, (b) the actual fields obtained by summing up 
modes. 
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Fig. 5.6 Plots of Poynting vectors of the actual fields of a 5X waveguide for 
various incident angles, (a) 0 = 0°, (b) 0 = 20°, (c) 0 = 40°. 





(a) RAY PICTURE 
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5.4 Summary 

This chapter demonstrated that the ray-optical description of the field in an open-ended 
waveguide excited by a plane wave is valid if the waveguide separation is large compared to the 
wavelength. As the waveguide separation gets smaller, the actual energy flow inside the guide 
begins to blur after propagating in a beamlike manner some distance into the guide. The simple 
ray-optical picture begins to break down at this point. Therefore, the GO-based SBR should 
be valid for deep cavities when the aperture is large compared to the wavelength. When the 
cavity is not too deep, SBR should also be valid for small apertures. For the small aperture, 
deep cavity case, the simple ray description must be modified to account for the beam blurring. 
This would possibly correct the discrepancy of the null predicted at 0 = 10° by SBR not found 
by modal anaysis for the <))(j>-polarization in Fig. 3.11(a). The solution to this problem will 
enable the extension of the ray-bouncing technique to a wider class of problems. 
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CHAPTER 6 

CONCLUSIONS AND EXTENSIONS 

In this thesis, the problem of the reduction of the radar cross section of open-ended 
cavities was studied. The issues investigated were reduction through lossy coating materials 
on the inner cavity wall and reduction through shaping of the cavity. A new method was 
presented to calculate the RCS of any arbitrarily shaped structure in order to study the shaping 
problem. The limitations of this method were also addressed. 

Chapter 2 studied modal attenuation in a multilayered coated waveguide. It was shown 
ihai by employing two layers of coating, it was possible to achieve an increase in both the 
magnitude of attenuation and the frequency band of effectiveness. It was also shown that the 
near axial incidence RCS of a terminated circular waveguide with a 4A. diameter aperture can be 
significantly reduced with a double layer of coating. As the angle of incidence increased, the 
effectiveness of the coating diminished. This is due to the fact that higher-order modes that are 
not highly attenuated by the given coating geometry are strongly excited at the wide angles of 
incidence. Further studies should be made in multilayered coatings to determine geometries 
that would result in a strong level of attenuation for higher-order modes in addition to the low- 
order modes. This would result in an increase in the range of effective RCS reduction not only 
in incidence angle, but in frequency as well, since more modes are excited in the guide as the 
frequency increases. 

The numerical method used in finding the roots of the characteristic equation breaks 
down when the coating thickness is very lossy and large in terms of wavelength. For complex 
arguments with large imaginary components, the Bessel and Neumann functions are related by 


J ^)(z) = iN^)(z) 
J'y(z) = iNVz) 


(6.1) 
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where i = and thus are no longer independent (shown in Appendix A). This causes the 
corresponding 4x4 matrices for that coating layer to become singular, and the matrix inversion 
in (2.9) is no longer possible. A reformulation of this problem using a different set of basis 
functions for the radial dependence of the modal fields might circumvent this problem. 

Chapter 3 presented a new method of computing the RCS of an arbitrary cavity, and 
Chapter 4 applied this method to study the effects of longitudinal bending on RCS reduction. 
Of course, experimental results to verify the SBR calculations would be highly desirable. For 
a complex structure such as the S-bend, SBR requires a large number of rays per wavelength 
for numerically convergent solutions. This results in a large computation time. The 
development of the SBR computer code and the results generated for this thesis would not have 
been possible without the availability of the NSF Cray-XMP supercomputer. A more 
sophisticated scheme in determining the far fields from the aperture field of the cavity might 
reduce the number of rays needed for convergent results and greatly reduce the computer time. 

Additionally, the SBR method could be applied to computing the RCS’s of external 
structures in addition to cavity structures. One such use would be to use SBR to compute the 
RCS of, for example, a realistic model of an entire aircraft instead of idealized models of 
different components of the aircraft. This has the added advantage of being able to take into 
account mutual interaction between the various parts of the aircraft. Another possible 
application of SBR would be for the analysis of beam waveguide problems [41]-[44]. 

Finally, Chapter 5 compared the ray and modal descriptions for the fields in a parallel- 
plate waveguide. For large waveguide separations, the fields do indeed exhibit a beam 
behavior after travelling a very large distance into the guide. For small separations, the beam 
picture begins to blur a short distance into the guide. To extend the range of validity of the 
SBR method, the simple ray picture must be modified to account for the beam blurring. One 
improvement would be launching of edge diffracted rays. This has the difficulty that, in 
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general, one ray gives rise to an infinite amount of diffracted rays. Another possible 
improvement lies in finding a proper way to account for the propagation of a finite-width 
collimated beam. 
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APPENDIX A 

BESSEL FUNCTIONS FOR LARGE ARGUMENTS 
WITH LARGE IMAGINARY PARTS 


For very lossy coating in an overmoded waveguide, kpr is complex with 

IkpTl »1 
Im(kpr) < 0 
llm(kpr)l » 1 

From Abramowitz and Stegun [45], (9.1.2) and (9.1.3): 

For Izl — > oo: J ^(z) = a /— cos (z - ^071 - 77 c) 

\ 7C 2 4 

N\)(z)= -v/^sin (z- 7 t) 7 t - 77 c). 

\ 7C 2 4 

Letz' = Z-— \) 7 t - — 7C. 

2 4 

If ”0 is fixed and finite, then 

Iz' 1 = |z --'UTt --Til as lzl-»oo. 

‘ ' 2 4 


Express z' as z' = x + iy where x and y are purely real. 
From Churchill [46], pg 57: 


. , . ,ey+e-y, ^ , 

sin z = sin X ( ) + 1 cos x ( 

2 

, ,ey+e-y. . . , 

cos z = cos X ( ) - 1 sin X ( 


ey-e-y 

2 

ey-e-y 


) 

). 


(A.1) 


(A.2) 


(A.3) 


(A.4) 


If y < 0 and lyl »1, then e*y »> ey, and we have 
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sin = 2 ^ ^ 

cos ^'=2 ®'^ X + i sin x). 


Clearly, cos z' = i sin z'. Therefore, 


J xXkpr) = iNx)(kpr) 
J'^/kpr)= iNVkpr) 


(A.5) 


(A.6) 


when (A.l) holds. 
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